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We propose a model for isotropization and corresponding 
thermalization in a nucleon system created in the collision of 
two nuclei. The model is based on the assumption: during the 
fireball evolution, two-particle elastic and inelastic collisions give 
rise to the randomization of the nucleon-momentum transfer 
which is driven by a proper distribution. As a first approximation, 
we assume a homogeneous distribution where the values of the 
momentum transfer is bounded from above. These features have 
been shown to result in a smearing of the particle momenta 
about their initial values and, as a consequence, in their partial 
isotropization and thermalization. The nonequilibrium single- 
particle distribution function and single-particle spectrum which 
carry a memory about initial state of nuclei have been obtained. 



1 Introduction 

The problem of isotropization and thermalization in 
the course of collisions between heavy relativistic ions 
attracts much attention, because, while describing 
experimental data, the application of thermodynamic 
models is one of the basic phenomenological approaches. 
Recently [1], this issue was examined for quark-gluon 
plasma produced as a result of ultra-relativistic A-\- A 
coUisions in experiments at CERN (Geneva) and BNL 
(Upton) [2]. 

Nowadays, there exist a few interesting models 
for the explanation of the thermalization phenomenon 
in parton systems. Among them, there is a model, 
which considers, in the framework of quantum 
chromodynamics, the influence of external color flelds 
on vacuum with the following creation of thermalized 
system of hadrons [3]. There is also a model 
of thermalization in heavy ion collisions, which is 
considered as a consequence of the Hawking-Unruh 
effect [4]. 

Such an enhanced attention to those phenomena 
arose, because, along with other factors, the assumption 
about a local thermodynamic equilibrium is successfully 
applied in various domains of high-energy physics. 



For instance, the transverse spectra of hadrons look 
thermalized not only if heavy ions collide, when the 
creation of many-particle statistical systems is adopted 
as an indisputable fact [5], but also at an e+e" 
annihilation (see work [6]), when the issue concerning 
the creation of many-particle system remains open for 
discussion. 

In this work, we attempted to flnd an explanation 
for the thermalization phenomenon, starting from such 
basic concepts as the conservation laws of energy and 
momentum, and leaving details, which are characteristic 
of every specific process, aside for further investigation. 
We offer a model, to a large extend simplified, which 
we called the maximal isotropization model (MIM). This 
model belongs, to a certain degree, to transport ones; we 
consider the evolution of the system, but we parametrize 
this development by the number of coUisions of every 
particle in the system rather than by the time variable. 
The idea of maximal isotropization consists in that, we 
suppose that a nucleon momentum transfer after every 
three nucleon-nucleon (-hadron) coUisions becomes a 
random quantity driven by a proper distribution. 

Let us assume for the moment that elastic scattering 
is the main contribution to the two-nucleon collision 
amplitude. Owing to the features of collisions between 
heavy ions, the initial (before the first collision) momenta 
of particles in A^-particle system A and A^-particle 
system B are known exactly. More precisely speaking, 
the initial momentum of every nucleon in the nucleus A 
is fca = feo = (0, 0, fcoz), while the initial momentum 
of every nucleon in the nucleus B is ki, — — fep — 
(0,0, —koz)- The energy and momentum are conserved 
in every separate coUision of two particles 

Uj{ka)+Uj{kb) = Uj{p^)+Uj{p^), ka + kb ^ Pa+Pb, (1) 

where ka and fcf, are the initial momenta of this 
particle pair, while p^ and Pb are the corresponding 
final momenta. We assume that the particles are on the 
mass shell, so that uj{k) = \/ -I- k^ . We also adopt 
the system of units, where ti — c— 1. 
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So, after the first collision, we have only four 
equations for the determination of six unknown 
quantities, and p^. It means that two quantities, e. g., 
{Pa)x and {Ph)x, remain uncertain and can be considered 
as such which accept random values, which are driven 
by a scattering probability. After the third collision, 
every component of the particle momentum becomes 
completely random. If the initial momentum is fixed this 
means a full randomization of the momentum transfer 
after three serial collisions. So, if we follow nucleon 
sequential elastic scattering from the first collision to 
the last one we would see a full randomization of 
the momentum transfer after every three sequential 
scattering. 




N N 

Fig. 1. Example of inelastic collision of two nucleons with 
creation of two 7r-mesons, N + N^N + N + tt + tt. 

In the inelastic collisions, the nucleon momentum 
transfer undergoes even faster randomization. Indeed, 
let us consider, for example, the process of creation of 

pions in the nucleon-nucleon reaction: N + N ^ 
N + N + ■ TT. The process leads to randomization 
of [3(2 + n,r) — 4] degrees of freedom, or dp degrees 
of freedom per particle becomes random, where dp = 
3 — 4/(2 + TT-Tr) and we assume that all particles are on 
the mass-shell. For instance, if n,r = 2 (see the physical 
diagram of the coUision depicted in Fig. 1), then two 
components of the momentum of every particle after 
reaction becomes random, i.e. dp — 2 in this process. 
If ^ 1, the number of the random degrees of freedom 
per particle achieves its maximal value dp = 3. We see 
that in the inelastic collisions the randomization of the 
nucleon momentum transfer, attributed to one physical 
diagram, goes faster than in the elastic scattering. 

To estimate the effective number of collisions, M, we 
analyze all nucleon collisions (physical diagrams), iVcoii, 
and obtain a total sum of the random degrees of freedom, 
c?tot, gained by the particle during all reactions before 
the freeze-out, i.e. we need to know dtot = SiJi" 
Then, the effective number of collisions is determined as 
M = dtot/3. 

So, in the framework of the MIM we assume a 
randomization of the momentum transfer, p (see Fig. 2), 



after every effective collision and the number of effective 
collisions, M, is determined for a particular experiment. 
At the same time, we assume that the value of the 
momentum transfer is bounded from above by a maximal 
value Pmax- Here, we based on the experimental fact that 



max 




Fig. 2. Nucleon momentum transformation as a result of 
two-particle collision; fc is the initial momentum, p is the 
momentum transfer, fc-l-p is the momentum after reaction, 
Pmax is the maximally allowed value of the momentum 
transfer. 

in hadron-hadron and nucleus-nucleus collisions at high 
energies the big values of the momentum transfer are 
suppressed and the eikonal approximation is appropriate 
for description of two-particle scattering and multi- 
particle production. 

2 Particle Distribution in Momentum Space 

As a building block for more compHcate picture of 
relativistic nucleus-nucleus collisions we consider an 
idealized many-particle system where every particle 
experiences the identical number of effective collisions, 
M, before system decay. Indeed, finite number of 
collisions is a common feature of the system with 
finite life time. Consider successive variations of the 
momentum of the n-th nucleon from nucleus A, which 
moves along the collision axis from left to right. Every 
TO-th coUision induces the momentum transfer, plT'K for 
the n-th nucleon so that, after M collisions, the nucleon 
acquires the momentum fc„: 

fco ^ ko+pl^^ ko+pH^+p'^^ > ko+Pn = kn , 

where P„ = J2m=i ) see Fig. 3. 




Fig. .3. Transformation of the initial nucleon momentum, 
fco, as a result of M collisions; P = II]^=i p''"'' is the total 
momentum transfer after M effective collisions. 
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To be more illustrative, the set of random variables 
p^™'' can be presented in the form of a table, where every 
row appears after every of the N particles from system 
A (the left side of the table) and B (the right side of the 
table) has collided: 





1 




N 


iV + 1 




2N 


1 

2 


fco 

p'? 




fco 


-fco 

„(i) 
Pn+1 

(2) 

Pn+1 




-fco 

„(i) 

P2N 

(2) 
P2N 


M 








PjV+l 




(M) 
P2N 



The final momentum fci of the first particle can be 
determined by summing up quantities in the first column 
of the table, the final momentum fc2 of the second 
particle by summing up quantities in the second column, 
and so on. 

The main goal of this work is to determine f2N, the 
density distribution function in the momentum space, 
which describes 2A'' nucleons after M collisions per 
particle. Wc write the total energy and momentum of 
the nucleons after freeze-out: 



2JV 



2JV 



and Pt, 



(2) 



n=l 



where e„ = uj{kn) = \J vn? + fcf^. All consideration 
is carried out in the c.m.s. of two identical coUiding 
nuclei, hence the initial total momentum of the system 
A + i? is equal to zero. As we discuss further, the 
total momentum of the system of 2N nucleons, Plot, 
can slightly differ from zero after freeze-out because of 
creation of secondary particles. Let us write down the 
density distribution function in the form 



f2N = Cf2N , 



(3) 



where C is the normalization constant. The 
unnormalized distribution function f2N can be defined in 
a two-fold way: first we follow all collisions of particular 
nuclcon by integration with respect to all nucleon 
random momentum transfer, second we fix the total 
energy and total momentum of 2A'^-nucleon system after 
freeze-out in a micro-canonical-Hke way. Then, it reads 

/2iv(-Etotj -Ptot; fel, • • • , fe2iv) = 



-I 



' dPi 

V 



2N 



V 



N 

n 



M 



6Hkn-ko-Y,P\ 



,(m) 



2N 

n 

n=N+l 



x5 Eu 



' M 

fc„ + fco - ^ P^n 
\ m=l 



2N 
n=l 



2N 



n=l 



(4) 



where V is the volume of the system in the coordinate 
space and the clement of volume in the momentum 
transfer space accessible for the n-th particle in the series 
of M coUisions reads 

dP„ ^ J (p(l)) d%'^ ..... J (p(f ') d'p'n''^ , (5) 

where the distribution of the momentum transfer 
is characterized by the presence of the form-factor 
J(p). For the sake of simpHcity we assume that this 
distribution is homogeneous one and the available values 
of the transferred momentum are restricted from above. 
So, as a first approximation we adopt that the form- 
factor has the following form: 

-''(P) = T7 6'(Pmax-|Pa;|)6'(pmax-|Pj,|)6'(Pmax-b2|), (6) 
Vp 

where Vp = 8p^^^ and J d^pJ{p) = 1. 

If the value of the maximally allowed transferred 
momentum scales as Pmax ~ 1 GeV/c then, 
our assumption actually means that nucleon-nucleon 
(hadron-hadron) reactions cannot go in the spatial 
volume which size is less than 0.2 fm or the nucleons 
cannot be closer than the distance 0.2 fm (compare this 
with the models where nucleon-nucleon potential has a 
hard core). 

Because we fix the total momentum of the nucleon 
system after freeze-out by the presence of the (5-function, 

(^Piot — fc?),^ , we introduce the multiplier Vp 

just to provide the dimension of the density distribution 
in the same way as of common use in statistical 
mechanics. 

Note, because of creation of secondary particles the 
energy of the nucleon system after freeze-out, .Btot, is 
in the following correspondence with the initial energy, 
Eq, which is fixed as the energy of net nucleons before 
collision in the c.m.s. of colliding nuclei 



Eq = Etot + Et 



sp ) 



(7) 



where E'sp is the energy of the secondary particles. 
And in a similar way the total momentum of the 
nucleon system after freeze-out, Ptot, is in the following 
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correspondence with the initial momentum, Pq, of the 
nucleons 



Po 



(8) 



where Pgp is the momentum of the secondary particles. 
If we investigate the coUisions of identical nuclei then, 
Po = 0. However, the value of the momentum of the 
nucleon system after freeze-out, Ptot, can be finite and 
is fixed by equation Ptot + Psp = 0. 

We normalize the density distribution function in 
such a way that simultaneously determines the density 
of states in the system: 



2N 



dki .. . dk2N /2w(-Etot, Ptot; fel, . . . , k2N) , (9) 



where the element of single-particle phase space looks 
like (in units of h): dkn = V ^2nyi ■ Hence, the 
normalization constant C is equal to 



1 



(10) 



Making allowance for notation ([5]) , the unnormalized 
density distribution ([4]) can be written down in the form 



f2N{Etot, Ptot; fei, 

/ 27V \ 



, k2N) 



2N 



n 

n=l 



V{2nf 



N 



■,ot ^ ^ k, 
I \ n=\ 

M „ 

,-io„-(fc„-feo) / 



use of this function, the unnormalized distribution 
density, /2Ar, can be rewritten in the unified form 



/2Af(-Etot, Ptot; fci, 

/ 2N 



, fc27v) 



27V 



N 

n 

n=l 
2N 



n=l 

1 f d^a. 



']]_ -ia„-(fc„-feo)+Aflng(o„) 



V J (27r)3 



n 



1 



'ib„-{k„+ko)+M In g{b„) 



(13) 



V J (27r)3 

To obtain the partition function, we should make the 
Laplace transformation of the density-of-states function 
^ with respect to the variable i?tot (it is evident that 
the function Sl2Ar depends on the external parameters 
E'tot and Ptot): 

oo 

^2w(/3,Ptot) = J rf£^tote-'^^*°'r!2Ar(£;tot,Ptot) = 

-Emin 



dki . . . dk2N /2Af(-Etot, Ptot; fci, ■ • ■ , k2N)- 



(14) 



This procedure, after the order of integration has been 
changed, brings about the Laplace transformation of the 
unnormalized density distribution function f2N' 



F2Ar(/3,Ptot;fel,---,fc 



27V) 



27V 

X 

n=7V+l 



n 



V{2Tr) 



d^P: 



M 



(m) 



oo 

J di;tote-^'''°727v(£^tot,-Ptot;fei,...,fc27v), (15) 



(11) Taking Eq. ifTS)) into account, we obtain 



27V 



Here, we presented (5-functions, which correspond to the _ 
conservation of momentum in a series of two-particle F2Ar(/3, Ptot; ^ij • ■ • j fc27v) = <5"^ I Ptot ~ ^ k 
collisions, in terms of the Fourier integrals over the V "=1 

variables a„ and 6„. 

Let us define an auxiliary function 



5(a) 



d^P ia- p A sin (a,;Pmax) 

=11- 



aiP-n 



(12) 



7V 

n 

n=l 
27V 



1 f d^a. 



V J (27r)3 



n_ g-ia„-(fc„-fco) + Mlng(o„) 



n 



'ibn-(kn+ka)+M In g(b„) 



V J (27r)3 

the Lapla 

(determined by ([6])) has been taken into account. Making consider the possible dependence of Pmax on Ei 



(16) 



,i=7V+l'- 

where the finiteness of the allowed momentum space While doing the Laplace transformation, we did not 
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Therefore, in what follows, the value of Pmax will be 
taken constant, thus playing the role of a parameter. 

So, the partition function lfT4|) . which is analogue 
of the partition function in the canonical ensemble (in 
what follows, for simplicity we name it as the canonical 
ensemble), can be written now in the following way 

Z2Ar(/3, Ptot) = Jdki.. .dk2N^2N{P,PtoU ki, . . .,k2N)- 

(17) 

Following this line one can determine the distribution 
function in the canonical ensemble 



F2iv (/3,Ptot;fci 



,1«2n)'- 



F2Ar(/3, Ptot;fcl, • ■ • ,fc 



2N) 



As a matter of fact, the mean values of the physical 
quantities which are evaluated^with the help of the 
distribution function, f2N — f2N/^2N from are 
obtained when the total energy of nucleons after 
freeze-out, iStot, is fixed. Meanwhile, the mean values 
evaluated with making use of the distribution function, 
F2iv(/3, Ptot; fci, ■ • ■ , k2N), will be at fixed parameter f3. 

3 Partition Function for a Large Number of 
Collisions 

To calculate integrals over the variables a„ and 6„ on 
the r.h.s. of Eq. (fTBl) in the case where the number of 
collisions for every particle is large, M ^ 1, the saddle- 
point method can be appHed. It implies the use of the 
following approximation: 



(18) 



This expression is to be substituted to iflGl) . The 
analogous operation is fulfilled in relation to the 
integration variables 6„. Using the Poisson integral we 
make integration in l|16|) with respect to the variables 
a„ and b„. Then, the partition function ^T7\ for the 
canonical ensemble, can be written in the following way 



^2iv(/3,Ptot) 



no 

n=l L 
2JV r 

n © 



(27r)3 , 

3/2 r (fkn 



(fxe 



-iPtofX 



(19) 



n=N+l 



(27r)3 

3/2 r d^kn 



i(fe„-fco)' 



(2^) 



where the (5-function, which refiects the conservation law 
of total momentum of the nucleon system after freeze- 
out, is presented in the form of Fourier integral over the 
variable x, and the parameter a is defined as: 



(20) 



2Mp 



2 

max 



Note, if instead of the Cartesian coordinate system one 
makes evaluations in the spherical coordinate system in 
the space of momentum transfer, then instead of the 
form-factor ^ and the quantity a defined in (|20|) one 
should use expressions: J{p) = Oij)^;^^ — \p\)/Vp, where 
Vp = 47rp^a^/3 and a = 5/{2Mpl^^^), respectively. 
We define auxiliary single-particle functions 



Za{b)iP,x) = (^^) 



3/2 



d^k 
(2^ 



xexp[-/3u!{k) + ik-x-a{k^ko)'^] (21) 

for the systems (nuclei) A and B, respectively. With their 
help, the partition function Z2n{(3, Ptot) from (fT9|) can 
be rewritten as follows: 

Z2N = t\ L32.e-»Ptof^+JVl„.„(/3,x)+iVln.,(/5,x)^ (22) 

(27r)3 J 

This expression will be basic for further research in this 
work. 

3.1. Collision of two identical nuclei 

In the case of collision between two identical nuclei, all 
nucleons in the one nucleus have the initial momentum 
ko and all nucleons in the other nucleus have the initial 
momentum — fco- The calculation of partition function 
(|22l) will be carried out making use of the saddle-point 
method at N ^ 1. We expand the functions In Za{P,x) 
and hiz{,{f3,x) (see l(2T|) ) into series with respect to the 
variable x about the point a; = up to the second order 
inclusive: 

^, i OZa(P,X) 

Xi 



In Za{f3,x)Ki\nZa{f3,0)+J2 
d^za{(3,x) 



Za{l3,x) dXi 



x=0 



2 ^ 



1 



dxjdxi 



Za{P,x) 
1 dZa{f3,x) dZa{P,x) 



x=0 



z'^{(3,x) dxj dxi 
We define two functions (see jSI])) 

Za{f3) = Za{,P,Q) , Zb(/3) =Z6(/3,0) 



(23) 



(24) 
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Then, regarding Za{/3) and Zb(/3) as the single-particle 
partition functions, the terms in brackets on the r.h.s. 
of l(23|) can be taken as the statistical averages of ki, 
kikj, and so on. Indeed, we designate such quantities by 
angle parentheses (. . .)^ and (. . \, where the subscripts 
a and b mean that averaging is carried out making use 
of the single-particle partition function Za{P) or Z},{P), 
respectively: 



(■ • ■)a{b) 



_ {a/n) 



3/2 



d^k 



(...)e 



Thus, we can rewrite expression ([23]) in terms of the 
averaged momentum components: 

3 

\nZa(j3.,x) W InZaiP) +i^{ki)aXi - 



1 ^ 

— 2 (^{kikj)a — {ki)a{kj)a 



(26) 



It is evident that the analogous expression can be written 
down for particles from system B as well. 

Now, we introduce the correlation function of 
momenta, 



c: 



(r) _ 



{kikj)r {ki)r{kj)r 



(27) 



where r = a or 5, and write expression l(26|) for system 
A and the analogous one for system B: 

1 ^ 



2 ^ 



It allows us to rewrite partition function ([221) as follows: 



Z2^(/3)«zf(/3)z,^(/3)^ 



X exp 
where 



(fa 



-iPtot ■ x+i2N{k) ■ x-N ^ djXtJ 



1 /"^(a) 

2 



a 



(28) 



(29) 



Integration of the Poisson integral on the r.h.s. of 
Eq. ([28]) brings about a simpler result 



Z2Nm « 4'(/3)4'(/3) 



1 



4ttN 



3/2 



(det c) 



1/2 



X exp 



(30) 



where = -Ptot/2iV. 

From expression ([30p . one can determine the 
"two-source" single-particle partition function, which 
simultaneously concerns both systems A and B: 



z{(3)=Za{P) Zbif3) 



Attn 



3/2N 



I 



l/N 



(det c) 



1/2 



(25) X exp 



(Pt„t - (fc)) • • (p,,t - (fc)) 



so that we can write down that 



Z2n{P) = z^'iP) 



(31) 



(32) 



We would like to emphasize the validity of the 
following limit value: 



lim 



1 



Attn 



3/2N 



= 1. 



For instance, at iV = 10, 100, 200, and 1000, the 
corresponding estimations are (l/(47riV))'^/^^ — 0.48, 
0.90, 0.94, and 0.99. As a result, in the case Ptot = 
and provided that the particle number N is rather large, 
it follows from Eq. ([31]) that 

lim z(/3) = z4/3)zb(/3)e-<'=>-^"-<'=>. 

3.2. Calculation of correlation matrix Cij 

Consider now the calculation of the elements of 
correlation matrix for the collision of two identical nuclei 
in the center-of-mass frame, i.e. we adopt the initial 
momentum fco — (0,0, fcoz). 

We would like to recall the definition of the 
correlation matrix. 



(a) 



c: 



(6) 



{kikj)a + {kikj)ii — {ki)a{kj)a ~ {ki)b{kj)b 



(33) 



We note, the average of transverse components of the 
momentum is equal to zero, {kx)a{b) ~ i^y)a{b) ^ 
From definition ([25]) it is evident that the effective 
average momentum ([29]) is 



{k)^^(^{k)a + {k)b 

while 

{kxky)a(b) — {^xkz) a{b) 



0,0,i((fe)„ + (fc,)f,) 



{kyk 



0, 



(34) 
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i.e. all the non-diagonal elements of the correlation 
matrix are equal to zero. As a result, the matrices C^°'^ 
and are diagonal; hence, the matrix inverse to C is 
diagonal as well, namely. 



/ 



V 



'■^11 +^11 






^22 "'"^22 




^33 ~^^33 



Consider the two-source single-particle partition 
function ((3T|) in the case where the total momentum of 
the nucleon sub-system after freeze-out equals to zero, 
PjQj — Ptot/2A^ = 0. Bearing the structure of the 
correlation matrix in mind and taking the value of the 
average moment (k) into account, partition function l(3T|) 
can be rewritten in the form 



z(/3) = z4/3) zb(/3) 



/ 1 



nil (c^^ 



-a 



1/2 



3/2N 



1/N 



(35) 



On the other hand, the z-component of the average 
momentum is equal to 





r (Pk 




1 {2ny 




r (fk 


2zb{l3) J 


(27r)3 



-kz exp \—[3u}{k) — a{k — feo)^] 

kz exp[-f3u;{k)-a{k+ko)^]. (36) 



We notice that the change of the integration variable 
kz —kz in integral IpT]) at a; = makes it evident 
that Za{(3) — Zb{(3). Using this equality, we obtain, from 
Eq. dMl), that 



{kz}- I ^kzC-^-^"^ 



Here, the integrand is no more than a product of an 
even and an odd functions. Owing to the symmetry 
of integration limits, the result of integration vanishes. 
At last, we have {kz) = \ ((^2)0 + {kz)}j) — 0) i-^- 



{kz) = — {kz)},- This result was expectedly obtained 
for identical nuclei when Ptot equals zero. 

Taking into account that Za{f3) = Zb{P) and 
the matrices (7^°^ and C''-''-' are diagonal, we draw a 
conclusion that those matrices coincide, i.e. 



g(a) ^ g(b) ^ 



(37) 



As a consequence, we can simplify Eq. l(35|) and write 
down the ultimate expression for the single-particle 
partition function: 



lim z{p) = Za{f3) zb{(3) 

N—*oo 



(?)V ^"^P [-Mka)-a{ka~ko)'] 



d^kb 

(2^ 



(2^) 

exp [-(3uj{kb) - a{kb + ko)^] 



(38) 



where the approximation 



l/N 



1 



AttN 



3/2N 



(det C) 



1/2 



fair at large enough N, was taken into account. 

Hence, on the basis of the proposed approach, 
the following nonequilibrium "two-source" distribution 
function has been obtained: 



/(fc(j, kb 



(«A)^ 

' Zaif3)zbi/3) 



exp [-Puj{ka)-a{ka - fco)^] 



X exp [-(3LL!{kb) - a{kb + fco)^ 
where in accordance with (|24l) 

d^k 



Za{b)iP) ^ (^) 



3/2 



(2^) 



-,3w(fc)-a(feq:fco) 



(39) 



(40) 



One can see that two-source distribution function ([39)) 
demonstrates features associated with the central Hmit 
theorem, which manifest themselves in the availability of 
two Gaussians, symmetrically located in the momentum 
space (remind, the analysis is carried out in the c.m.s. of 
two colliding nuclei) . Actually, the latter is the reason for 
the name "two-source". The proposed approach allowed 
not only the expected general result to be obtained, 
but also the expression for the dispersion and the 
distribution centers to be deduced, which can be checked 
up experimentally in nucleus-nucleus collisions. 
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In particular, we note that the increase of the 
number of collisions M and the accessible volume in the 
momentum transfer space, which is determined by the 
parameter pmax (see l(20|) ). gives rise to the smearing of 
this effect and the transformation of the distribution into 
the thermal one. 

3.3. Single-particle spectrum 

Consider the spectrum of particles, SN/d^p, after 
freeze-out. If every nucleon from the large system, 
according to the proposed model, takes part in M 
collisions, the spectrum can be obtained by averaging 
the quantity (as well as an arbitrary function D{ka, kb) 
of random variables) 



D{p, ka, kb) = NS^ip - ka) + NS^ip - kb) 



(41) 



over the values of fca and kb with making use of two- 
source single-particle distribution function ((39l) . i.e. 



(27r)3 (27r)3 



D{p, ka,kb) f{ka,kb) . 



(42) 



Then, after averaging the random quantity JS]), the 
spectrum, d?N/d?p — (D), looks like 



(PN 



N 



3/2 



e~/3-(p) X 



1 



g-a(p-fco)^ _j_ 



1 



-Q(p-l-fco) 



(43) 



MP) ^fc(/3) 
with a defined in |[20]). The spectrum |[43]) is 
obtained as a sum of two terms which represent 
the two-source structure: the first contribution oc 
exp [—f3uj{p) — a{p — fco)^]/^a is associated with nucleus 
"A", and the second oc exp [—(}uj{p) — a{p + feo)^]/zh is 
associated with nucleus "B". 

If we consider the initial momentum fco — (0, 0, /coz), 
the spectrum possesses two maxima near the two values 
of the longitudinal momenta: pz — fcoz and pz — —koz. 
Let the notation p^ stands for p^ +pl; then, the particle 
spectrum in the momentum space looks like 



(Pp 



zoiP) 



g-/3a)(p)-Qp^ X 



X 

where 



g-Q(pj-fcoj) _|_ g-a(P2+fco2) 



Vvr/ J [2n)'^ 



(44) 



(45) 



3.4. Effective temperature at kg = 

It is worth noting that in kinetic approach, where 
the distribution function depends on time, f{t,k), the 
temperature can be determined just conventionally. The 
latter problem is analogous to a definition of instant 
velocity of the object, which moves with acceleration. 
Actually, in our approach we parametrize the time axis 
by the number of particle collisions M = {v)t, where 
{v) is the mean frequency of collisions and t is the time 
interval. In this sense the obtained distribution function 
(|39l) with respect to variable M can be regarded as 
a nonequilibrium one (M oc time) and the parameter 
P = 1/T cannot be regarded as an inverse temperature. 
Just in the asymptotic regime when t — » 00 or on 
our language M ^ 00 the distribution function (|39| 
transforms to the Boltzmann one and the parameter T 
becomes the Boltzmann temperature. Note, for the sake 
of simplicity we do not consider special statistics so far. 

Model, which describes a subsystem of particles with 
a zero initial momentum, fco = 0, can be applied as a first 
approximation to the secondary particles, for instance 
pions, which are created in relativistic nucleus- nucleus 
collision. One can regard the momentum of the particle 
just after creation as a first momentum transfer, which 
is added to zero initial momentum. On the other hand, a 
particular example discussing in the present paragraph 
can be regarded as a mathematical limit when fco ^ 0. 

First, consider particles, which are created near 
threshold. In this case a nonrelativistic behavior is 
relevant to single-particle energy uj{p) w 171+ p^ /2m, 
where m is the particle's mass. The subsystem under 
consideration consists of N particles: for condition fco = 
the distribution function II39I) reduces to 



(46) 



z{(3) 



where 



Evidently, for this initial condition the two-source 
distribution (|39|) reduces to the "single-source" one. 

An estimation for the effective temperature can be 
obtained in the following intuitive way: 



g 2mT "'^ 



(27r)3 

where 
11 1 

eft J J coll 



mT, 



cff 



27r 



3/2 



(48) 



(49) 
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and the collision "temperature" is defined as 

Mpl 



T, 



coll 



52 

max 



3m 



(50) 



For the thermal wave length determined as A = 
yJiir/mT and, respectively, Aoff = y/2ji/rnT^, the 
spectrum can be written down in the form 



(Pp (27r) 



^ exp 



2mT, 



off 



(51) 



We conclude immediately: the restriction of the 
accessible volume in the space of transferred momentum 
and the finiteness of the collision number turn out to 
effectively reduce the temperature, Teff < T. Really, 
according to Eq. l(49|) . 



T 



1 + r/TeoU 



(52) 



Thus, the increase of the collision number M is 
accompanied by the growth of the effective temperature 
Tcff up to its limit value T. It is the reason of why, 
when the parameter M is large enough, the quantity 
Toff in Eq. ijSTI) should be replace by T. On the other 
hand, if Pmax tends to infinity, spectrum ((5T|) acquires 
the standard form 



lim 



In the general case (T 



(53) 

m) , the relativistic 
f has to be 



dispersion relationship (j-'(p) — ^Jrn? 
used. Let us make the inverse Laplace transformation 
of the many-particle partition function: 



f^Af(-Etot) 



1 

27ri 



C+ZOO 



(54) 



In order to calculate the integral on the r.h.s., we use 
the saddle-point method. For this purpose, the last 
expression is rewritten in the form 



The condition for the extremum to take place along the 
real axis (the variable (3 is real) looks like 



If 1 dz{p) 

N z{l3) dp 



(56) 



The solution of this equation with respect to /3 gives the 
value of the parameter T = 1/(3, which corresponds to 
the average energy value per particle, namely, Etot/N. 
Here, we may use the explicit form of the single-particle 
partition function l(47|) to obtain 



1 _{aM^ 
N ~ z(/3) 



(27r)3 



uj{k) e 



(57) 



We obtained a transcendental equation with respect to 
the parameter (3. The value od this parameter, which 
is the root of equation, differs from the value, which is 
determined from the similar equation for the ideal gas 



N Zid(p) J (27r)'^ 



where 



(58) 



(2^ 



-f3ui{k) 



(59) 



is the single-particle partition function for the ideal gas. 

Similarly to what was done in the non-relativistic 
case, one can introduce the idea of effective temperature 
Tcs = l//?off in the following way: 



1 

27ri 



c+ioo 



d(3e 



F{0) 



(55) 



1 dz{(3) 
8(3 



1 



dzidiPcs) 



where Eq. I|32l) was used and we introduce notation 
F{(3) = (3Etot + N\ogz{(3). Now, the minimum of the 
function F{(3) along the imaginary axis and its maximum 
along the real axis of the variable /? are to be determined. 



where in the non-relativistic case the equality is valid 
when /3eff = /3 + l/Tcoii- I.e., by definition the effective 
temperature, is the temperature of ideal gas, at which 



9 



D.V. ANCHISHKIN, S.N. YEZHOV 




50 100 150 200 
T (MeV) 



% 500 




100 150 
T (MeV) 



200 



Fig. 4. The average single-particle energy after 5 (upper 
panel) and 10 (lower panel) consecutive effective scattering 
for Etot per particle equals 200, 400, and 800 MeV (solid 
curves from bottom to top); and the average single-particle 
energy for ideal gas of vr-mesons (dashed line). 

the average single-particle energy, which is equal to 
Etot/N, is the same: 



(a/7r) 



3/2 



z(T) 



V 



Zid(T'cff) J (27r) 



uj{k) e 



Lu{k) e 



(60) 



In our illustrative example we assume that the pion 
energy, which corresponds to the maximal momentum 
transfer, pmax, scales like the mean energy per particle: 
Etot/N — y/ -|- 3p^ax : where is the pion mass. 
In Fig. 4, the results of calculations of single-particle 
average energy of pions in accordance with Eqs. I|57p 



(solid curves) and (|58l) (dashed curve) are depicted as 
the function of the parameter T (temperature for the 
case of ideal gas). The horizontal straight line in this 
figure corresponds to a constant value of single-particle 
average energy, {uj{k)) — const. The intersections of this 
straight line with a curve determine the temperature, 
which corresponds to this value of the single-particle 
average energy. Therefore, fixing the average value of 
the energy per particle, we come to a conclusion that 
the effective temperature determined by the intersection 
with the dashed curve is lower than the value of the 
parameter T, which is determined by the intersection 
with the solid curve, i.e. Toff < T. 

Such a situation completely corresponds to the 
relationship between the effective temperature and 
parameter T in non-relativistic systems (|52l) . Similarly 
to the non-relativistic case, the effective temperature 
is a limit value for parameter T, which is reached 
after a sufficiently large number of collisions M occurs 
and/or provided that the accessible momentum space is 
unconfined, i.e. T Toff, if Af ^ c« and/or pmax —>■ oo. 

4 Discussion and Conclusions 

In this work, on the basis of the model proposed for 
the collision of two relativistic many-particle systems 
(for example, nuclei), the process of thermal equilibrium 
establishment has been analyzed. This process was 
parametrized by the number of effective collisions per 
nucleon: every effective collision brings about the total 
randomization of the nucleon momentum transfer. Note 
that in kinetic models, in contrast to our approach, 
the processes of reaching the thermal equilibrium are 
parametrized by time rather than the collision number. 

The main result of the work is the explicit expression 
(|39|) for the two-source nonequilibrium distribution 
function, which depends on the collision number and the 
effective volume of the single-particle space of allowed 
momentum transfer. This volume is determined by the 
parameter Pmax- There is one more essential feature of 
the distribution function (|39| : it carries memory of the 
state of nuclei before collision, i.e. the initial momentum 
of nucleons, feo, of one nucleus and the initial momentum 
of nucleons, — feo, of another nucleus. The immediate 
successful application of the two-source distribution ((39l) 
is evaluation of the single-particle spectrum l(43|) . The 
spectrum has two items, which can be attributed to the 
first and to the second colliding nuclei, respectively, and 
in accordance with this structure it can be named as a 
"two-source single-particle spectrum". 

In application to relativistic nucleus-nucleus 
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collisions we extend a statistical model where the many- 
particle partition function, which describes the many- 
particle system just after freeze-out, is defined as a 
multi-dimensional phase-space integral (see e.g., [8-11]). 
In this approach the single-inclusive particle spectrum, 
which is created in multiparticle production process, 
is attributed to the volume of 3(iV — 1) dimensional 
momentum space of the unobserved — 1 particles 

w(fci) oc |M(i^,fci,fc2, ■ • ■,kN)f 

where /ci, . . . , fc^r are the final momenta of all particles 
after freeze-out with fc° = w(fci), K is the final 
total 4-momentum of the system of these particles, 
and k2,---,kN are the fixed values of the momenta 
determined by the law of the mean when one extracts 
\M{K,ki,k2, . . . ,kis[)\^ from the integrand on the r.h.s. 
of (|6T1) . Note, |M|^ is the modulus squared matrix 
element of the multiparticle production process, which 
in a simplest version of the statistical model is adopted 
as a constant due to phase-space dominance. 

Actually, the matrix element of the multiparticle 
production, M{K,ki . . . ,kN), does not describe the 
interactions of the created particles (final-state 
interactions [10]) before decay of the system. At the same 
time, these interactions are responsible for a partial or 
a full thermalization of the system. 

In future, in the spirit of the present work, we plan to 
include in l(6T|) the rescattering of the secondary particles 
after their creation till the freeze-out in the same manner 
as it is proposed above within the maximal isotropization 
model. 

So, in the framework of the quasi-microcanonical 
concept (total energy is fixed) of the statistical model, 
l(6T|) . one cannot say how many rescattering every 
particle went through after creation and before freeze- 
out, moreover if there been any rescattering at all 
(imaging that the system of neutrinos is created). 
Actually, without answer to this question it is not 
possible to shift consideration to the canonical ensemble 
(temperature-fixed) [11]. That is why, with respect just 
to expression (|6T1) and without complementary kinetic 
arguments, which account for particle rescattering, one 
cannot say to what extent thermalization process takes 
place in the system before the particles are freezed out. 

At the same time, in our approach we considered 
successive collisions between particles during evolution 



of the system and obtained the nonequilibrium single- 
particle distribution function l(39|) . which shows the 
degree of thermaHzation in the system. That is, 
we considered the process of thermal equihbrium 
establishment as a function of the finite number of 
collisions per particle, M, before freeze-out. It has to 
be noted that the parametrization of a nonequilibrium 
process by the number of coUisions per particle is 
especially promising in the range of relativistic energies, 
because the coUision number is an invariant with respect 
to the change of the reference frame. 

We parametrize the time axis by the number of 
particle collisions M — (i')i, where {v) is the mean 
frequency of coUisions and t is the time interval. If one 
considers a dependence of the distribution function (|39l) 
on the variable M, exactly in this sense he can regard 
this distribution as nonequilibrium one (M oc time). 
Increase of M effectively mimics an approach to the 
Boltzmann distribution (for the sake of simplicity we 
do not consider special statistics in this issue). If the 
number of collisions, M, is fixed, it means that the 
process of thermalization stopped at the moment when 
the particles had experienced just M collisions what 
results in a partial thermalization of the system. So, the 
level of thermalization and isotropization depends on the 
number of effective collisions, M, which is determined, 
first, by the life time of the system or more precisely by 
the number of physical collisions of every particle, iVcoii, 
and, second, by the level of randomization of momentum 
transfer in every physical coUision. As we have discussed, 
a degree of randomization of the particle momentum 
transfer attributed to a single physical diagram (see 
Fig. 1) varies from one for elastic collisions to three for 
inelastic reactions, when a yield of secondary particles 
is large. Note, this concerns just the hadronic processes, 
and we keep aside so far all issues related to the quark 
gluon plasma. 

So, we understand the relation of the value of the 
collision number M to the degree of thermalization in 
the system. Let us discuss how the number of effective 
coUisions depends on coUision energy. For instance, we 
make a shift from AGS energies to SPS ones (we fix 
centrality and the type of nuclei). First, with increase of 
energy the number of physical hadron collisions, A'coii, 
decreases. Indeed, when one considers the sequential 
reactions of a nucleon in the framework of UrQMD 
[12, 13], he finds a less number of the physical diagrams 
(like that depicted in Fig. 1) at SPS energies than at 
AGS energies. Remind, the number of effective collisions 
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is determined as: 

^^=^E4^ (62) 

^=l 

where dp ^ is the number of random degrees of freedom 
gained by every particle during i-th coUision (see 
discussion after Fig. 1). Hence, with increase of collision 
energy the number of items in the sum (|62| decreases. 
On the other hand, the mean value of every item in the 
sum ((62l) increases because of increase of the yield of 
secondary particles. However, as we showed, the creation 
of secondary particles can increase the number, M, of 
effective coUisions not more than by the factor three. So, 
there are two mechanisms, which determine the number 
of effective collisions M and hence they determine the 
degree of thermalization of the nucleon system. At 
the same time, these two mechanisms compete with 
one another with increase of the collision energy. For 
instance, in case of their balance, the dependence of M 
on the collision energy is absent. It is necessary to add for 
completeness, that the parameter Pmax, which influences 
also on the degree of thermalization (see definition (|20|) ). 
evidently increases with increase of coUision energy. 

To give an explicit example we evaluate the effective 
temperature of the subsystem of particles for particular 
initial condition: initial momentum of particles equals 
zero, fcg = 0. This case may be regarded as a toy model, 
which considers the rescattering of the particles created 
in nucleus-nucleus collision (momentum of the particle 
just after creation is regarded as a first momentum 
transfer) or as a mathematical limit when fco ^ 0. 
Effective temperature. Toff, is defined as a temperature 
of the ideal gas of the same particles with the same mean 
energy per particle. It turns out that the temperature of 
the ideal gas under this condition is always lower than 
the inverse parameter T = 1/ (3 from the nonequilibrium 
distribution (|39l) when the number of collisions, Af , is 
finite. The temperatures biggin to coincide, Tcff T, 
when collision number is big enough, M ^ oo. 

We hope that the distribution function derived in the 
present paper will help to understand better the effect 
of the hadronic rescattering in non-central relativistic 
collisions of heavy ions on the azimuth anisotropy of the 
momentum spectra (eUiptic fiow) [14]. 



A comparison of the results obtained theoretically 
with experimental ones will be carried out elsewhere. 
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